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Abstract. We review theoretical results obtained recently in the framework of sta- 
tistical mechanics to study systems with long range forces. This fundamental and 
methodological study leads us to consider the different domains of applications in a 
trans-disciplinary perspective (astrophysics, nuclear physics, plasmas physics, metallic 
clusters, hydrodynamics,...) with a special emphasis on Bose-Einstein condensates. 



1 Introduction 

Properties of systems with long range interactions are to a large extent only 
poorly understood although they concern a wide range of problems in physics. 
Recently, the disclosure of new methodologies to approach the study of these 
systems has revealed its importance also in a trans-disciplinary perspective (as- 
trophysics, nuclear physics, plasmas physics, Bose-Einstein condensates, atomic 
clusters, hydrodynamics,...). The main challenge is represented by the construc- 
tion of a thermodynamic treatment of systems with long range forces and by 
the understanding of analogies and differences among the numerous domains of 
applications. 

Some promising results in this direction have been recently obtained in the 
attempt of combining tools developed in the framework of standard statistical 
mechanics with concepts and methods of dynamical systems. Particularly ar- 
duous, but very exciting, is the understanding of phase transitions (negative 
specific heat, non extensive thermodynamics,...) for such systems as well as all 
the aspects related to non-equilibrium phenomena and their description in terms 
of dynamical concepts (self-consistent chaos, slow relaxation, formation and role 
of structures,...). 

Finally this fundamental and methodological study should help us to de- 
tect the depth and the origin of the analogies found in the different domains 
mentioned above or on the contrary emphasize their specificities. In particular, 
we would like to put a special emphasis on Bose Einstein Condensation (BEC) 
which could be the main field of applications, since experiments and theoretical 
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ideas have reached an impressive quaUty in the last decade. In that domain, 
many inequivalences between ensembles have been reported and should be clar- 
ified. Moreover, long range interactions in BEC have opened very exciting new 
perspectives to consider BEC as a model for other systems. 

2 Why systems with long range interactions are 
important ? 

2.1 The problem of additivity 

The methods to describe a given system of N particles interacting via a gravi- 
tational potential in 1/r are dramatically dependent on the value N. If Newton 
showed the exact solution for N — 2, and one can expect to get a numerical 
solution in the range = 3 — 10^, the results are clearly out of reach for a 
larger number of particles. In addition, it is clear that the knowledge of the evo- 
lution of the different trajectories is completely useless, since it is well known 
that these systems are chaotic as soon as N is greater than two. Therefore, one 
needs to get a statistical analysis, in order to get insights in the thermodynamical 
properties of the system under study. 

However, such statistical study leads immediately to unexpected behaviors 
for physicists used to neutral gases, plasmas or atomic lattices. The underlying 
reason is directly related to the long range of the interaction, and more precisely 
to the non additivity of the system. 

To avoid misunderstandings, let us first clarify the definition of extensivity 
with respect to additivity. A thermodynamic variable, like the energy or the 
entropy, would be extensive, if it is proportional to the number of elements, once 
the intensive variables are kept constant. To be more precise, let us consider the 
following Hamiltonian 

'' = -i{p) ' 

where the spins Si — ±1 sit on a one-dimensional lattice labeled by i = 1, . . . ,N. 
Without the 1/N prefactor such a Hamiltonian would have an ill defined ther- 
modynamic limit. This is correctly restored by applying the Kac prescription 
within which the potential is rescaled by an appropriate volume dependent fac- 
tor, here proportional to N since the lattice is one-dimensional: such a Hamilto- 
nian is therefore extensive. Let us note in passing that this regularization is not 
always accepted. In cases with a kinetic term, such a regularization corresponds 
to a renormalization of the time scale. 

On the contrary, this Hamiltonian is not additive. Indeed, let us divide a 
system schematically pictured in Fig. 0, in two equal parts. In addition, one 
considers the particular case with all spins in the left part are equal to 1, whereas 
all spin in the right part are equal to -1. It is clear that the energy of the 
two different parts, will be Ei ^ E2 ~ —j^ (y) — However, if one 



Dynamics and thermodynamics of systems with long range interactions 



3 



computes the total energy of the system, one gets = — ■^(■y — "y) =0. It is 
therefore clear that such a system is not additive, since one cannot considers that 
E1+E2 = E, even approximately. The energy of the interface, usually neglected, 
is clearly of the order of the energies of the two different parts: the system is 
not additive. The underlying reason is that Hamiltonian (jl|) is long (strictly 
speaking infinite) range, since every spin interact with all others: moreover, as 
the interaction is not dependent on the distance between spins, this is a mean- 
field model. This example is further elaborated in ji). 
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Fig. 1. Schematic picture of a system separated in two equal parts. 



This non additivity has strong consequences in the construction of the canon- 
ical ensemble. Once the microcanonical ensemble has been defined, the usual 
construction of this ensemble is usually taught as follows. The probability that 
system 1 has an energy Ei and system 2 has an energy E2 is proportional to 
fii{Ei) i72(£'2) dEi, where the number of states of a system with a given en- 
ergy E, f2{E), is related to the entropy via the classical Boltzmann formula 
S{E) ~ InE (we omit the factor for the sake of simplicity). Using the ad- 
ditivity of the energy, and considering the case where system 1 is much smaller 
than system 2, one can expand the term 6*2 (-E — Ei), as shown in the following 
different steps 

ni{Ei) f22{E2) dEi = Q^iEi) n^iE - El) dE^ (2) 
= ni{Ei) e^^(^ ~ ^1) dEi (3) 

iS2iE)~Ei^ +...) 
= Qi{Ei)e oE\E dEi (4) 

oc ni{Ei) e^^^i dEi , (5) 

where /? — ^^|^- One ends up with the usual canonical distribution. It is clear 
that additivity is crucial to justify (|^), which means that non additive system 
will have a very peculiar behavior if there are in contact with a thermal reservoir. 
This is one of the topic discussed in this paper, and in numerous contributions 
in this book. 



4 



T. Dauxois et al. 



2.2 Definition of long range systems 

To define now systems with long range interactions, let us consider the potential 
energy for a given particle, situated in the center of a homogeneous sphere with 
radius R. We will omit at this stage the interactions of the particles situated in 
the neighborhood of the particles by forgetting particles located in the sphere of 
radius e, where e R as shown in Fig. The reason will be explained in the 
following subsection. 



Fig. 2. Schematic picture of a particle interacting with all particles located in a homo- 
geneous sphere of radius R, except the closest ones located in the sphere of radius e. 

If one considers that particles interact via a potential energy proportional to 
1 /r" , where a is the key-parameter defining the range of interaction, we obtain 
in the three dimensional space 



where p is the particle density. When increasing the radius R, the contribution 
due to the surface of the sphere, could be neglected when a > 3, but 

diverges if a < 3. In the latter case, surface efi'ects are important and therefore 
additivity is not fulfilled. 

If one generalizes this definition to long range systems in d dimensions, one 
easily shows that energy will not be additive if the potential energy behaves at 
long distance as 



Mean-field models as Hamiltonian (|l|) corresponds to the value a = 0, since the 
interaction does not depend on the distance. There are therefore not additive 
as shown in section J. Barre et al consider Q such a mean field model: 

the Blume-Emery-Griffith (BEG) model with infinite range interactions. The 
gravitational problem, which is at the origin of this study, and corresponds to 
a = 1 in three dimensions, clearly belongs to this category, but presents also 
additional difficulties. 







< 1 



(7) 
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2.3 Difficulties with the gravitational problem 

This problem is particularly tedious because, in addition to the non additiv- 
ity due to long range character, such a system needs a careful regularization 
at short distances to avoid collapse. To be more specific, let us consider the 
configurational partition function of a system of N particles 

Zu^ f d'^rt , (8) 

Jv 

where one notes U{ri ) the gravitational potential energy, /3 the inverse of the 
temperature and V the volume of the system. From the shape of the potential 




Fig. 3. The gravitational potential energy as a function of the distance r is represented 
by the solid curve, whereas the dotted one shows the regularized potential energy to 
avoid gravitational collapse. 

energy depicted in Fig. ||, one clearly see that Zu will diverge if all particles 
collapse towards the same point. This difficulty arises because the potential 
energy is not bounded from below as for a Lennard- Jones or a Morse potential. 
This effect is of course physically forbidden by the Pauli principle. However, to 
avoid the use of Quantum Mechanics, the usual trick is to introduce an add-hoc 
cut-off. The potential is therefore "regularized" by introducing the value — C, 
shown in Fig. |[ Thus, the inequality U{Ti ) > —C allows easily to find a finite 
upper bound for the configurational partition function Zjj < V e'"^, where V is 
the volume of the system. 

However, there is a third difficulty in the case of gravitational interaction: the 
system is open, i.e. without boundary, strictly speaking. In the microcanonical 
ensemble , the number of states 

^{E)= jWd^P^j \{dq^5{E-H{q,,pi)) (9) 

i i 

will diverge if the system is not confined. This divergence is actually not re- 
stricted to the gravitational interaction but would also occur if one considers a 



6 



T. Dauxois et al. 



perfect gas in an infinite volume. However, one considers of course always a gas 
in a finite domain, i.e. in a finite volume. This is not any more possible for the 
gravitational interaction where the system is clearly infinite. 

Despite these additional difficulties, the astrophysics community has obtained 
an impressive quantity of results in this domain. Thanu Padmanabhan de- 
scribes several remarkable features, both for isolated gravitating systems as well 
as for systems undergoing nonlinear clustering in an expanding background 
cosmology. The emphasis is on general results and he brings out the inter- 
relationships of this subject with topics in fluid mechanics, condensed matter 
and renormalization group theory. 

Similarly, Pierre-Henri Chavanis |^ presents how the structure and the orga- 
nization of stellar systems (globular clusters, elliptical galaxies,...) in astrophysics 
can be understood in terms of a statistical mechanics for a system of parti- 
cles in gravitational interaction. Finally, Eddie Cohen and laroslav Ispolatov [Q 
consider the related gravitational-like collapse of particles with an attractive 
1/r" potential. Using mean field continuous integral equation, they determine 
the saddle-point density profile that extremizes the entropy functional. For all 
0<a<d = 3, a critical energy is determined below which the entropy of the 
system exhibits a discontinuous jump. 

2.4 Applications to large systems 

However a growing scientific community has recently begun to tackle this prob- 
lem with different viewpoints. One of the fascinating aspects of this problem is 
that, in addition to gravitating systems, it concerns a large variety of systems 
that we would like to discuss briefly in the following section. 

Plasmas : Rarefied plasmas share many properties with coUisionless stellar 
systems, and in particular that the mean field of the system is more important 
than the fields of individual nearby particles. Here again, the Coulomb force is of 
long range character. However, there is a fundamental difference between plasmas 
and gravitation. Plasmas have both positive and negative charges, so that they 
are neutral on large scales and can form static homogeneous equilibria; on the 
contrary, gravitating systems can never form static homogeneous equilibria. This 
so-called Debye screening explains why many techniques of plasma physics can 
not be transferred immediately. Elskens [|| and Diego Del Castillo Negrete 
present some of their results in the framework of plasma physics. 

2D Hydrodynamics : Two-dimensional incompressible hydrodynamics is an- 
other important case where the interaction is long range. Indeed, the stream- 
function tp is related to the modulus of the vorticity oj, via the Poisson equation 
Atp = uj. Using the Green's function technique, one easily finds that the solution 
is 
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where G (t' — T*'') depends on D, but G (T*) ~ | hi T^|, when ^ 0. The kinetic 
energy being conserved by the Euler equation (dissipativeless), it is straightfor- 
ward to compute it on the domain D, with boundary dD, 

E = jj^^\{^i^f (11) 

ltdl iljWtp +- 1 (f'^ uj{~r')il}{~t) (12) 

dD 2 

= / / d^T^d^T*' ti;(7^')cj(T*) hi |~?^ — (13) 

47r J Jd 

since ^ = on dD. This emphasizes that one gets a logarithmic interaction. 
The analogy is even more clear if one approximates the vorticity field by point 
vortices oj(~r') = ~ located at ~f*i, with a given circulation r^. 

The energy of the system reads now 

E^\Y.r.rM-^^-^',\ ■ (14) 

The interaction among vortices has a logarithmic character, which corresponds 
to a = 0. 

Pierre-Henri Chavanis |^ studies carefully the analogy between the statis- 
tical of large-scale vortices in two-dimensional turbulence and self-gravitating 
systems. This analogy concerns not only the equilibrium states, i.e. the forma- 
tion of large-scale structures, but also the relaxation towards equilibrium and 
the statistics of fluctuations. Diego Del Castillo Negrete discusses also his 
results in the framework of hydrodynamics. 

Dipolar interactions Dielectrics and diamagnets in an external electric or 
magnetic field exhibit a shape dependent thermodynamic limit pl| . This is due 
to the marginal decay of the potential energy a =^ d ~ 3 for systems of dipoles. 
There is some approach to the solution of this problem only in zero field and in 
the absence of spontaneous ferromagnetisni jl2j. This is a border case for the 
long-range interactions, but it deserves a special attention. 

Fracture Let us examine analytical solutions for the plane stress and displace- 
ment fields around the tip of a slit-like plane crack in an ideal Hookean continuum 
solid. The classic approach to any linear elasticity problem of this sort involves 
the search for a suitable "stress function" that satisfies the so-called biharmonic 
equation V^(V^?/') — where is the Airy stress function, in accordance with 
appropriate boundary conditions. The deformation energy density is then de- 
fined as U oc ae where a is the fracture stress field around the tip, whereas e 
is the deformation field. Considering a crack-width a in a two-dimensional ma- 
terial and using the exact Muskhelishvili's solution jlOj, one otbains the elastic 
potential energy due to the crack 

^ ^go(l-^) ^ 
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where E is the Young modulus, a^o the stress field at infinity, u the Poisson 
coefficient and r the distance to the tip: the elasticity equation in the bulk of 
solids leads therefore, again, to a border case for the long-range interactions 
since U ~ in two dimensions. It appears that, despite of its engineering 
applications, the dynamics of this non conservative system has been very little 
studied, presumably because of its long range character. In addition, in such a 
two dimensional material, the presence of several fractures could exhibit very 
interesting screening effects. 



Table 1. Table listing different applications where systems are governed by long range 
interactions. Large systems where the interactions is truly long range and small systems 
where the range of the interactions is of the order of the size of the system are separated. 



Interactions 


a 


a/d 


Comments 


Large systems 


Gravity 


1 


1/3 


Long range 


Coulomb 


1 


1/3 


Long range with Debye screening 


Dipole 


3 


1 


Limiting value 


2D Hydrodynamics 








Logarithmic interactions 


Fracture 


1 


1 


Stress field around the tip 



Small systems 



atomic and molecular clusters 
nuclei 

BE Condensation 



2.5 Applications to small systems 

In addition to large systems where the interactions is truly long range, one should 
consider small systems where the range of the interactions is of the order of the 
size of the system. In these cases, the system would not be additive either and 
many similarities will be enc;ountered. Phase transitions are universal properties 
of interacting matter which have been widely studied in the thermodynamic 
limit of infinite systems. However, in many physical situations this limit is not 
accessible and phase transitions should be considered from a more general point 
of view. This is for example the case for some microscopic or mesoscopic systems: 
atomic clusters can melt, small drops of quantum fluids may undergo a Bose- 
Einstein condensation or a super-fluid phase transition, dense hadronic matter is 
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predicted to merge in a quark and gluon plasma phase while nuclei are expected 
to exhibit a liquid-gas phase transition. For all these systems the experimental 
issue is how to characterize a phase transition in a finite system. 

Philippe Chomaz and Francesca Gulminelli discusses results from nu- 
clear physics as well as from clusters physics. In particular, they introduce a 
possible definition of first order phase transitions in finite systems based on 
topology anomalies iQ of the event distribution in the space of observations. 
This generalizes the definitions based on the curvature anomalies of thermody- 
namical potentials and provides a natural definition of order parameters. The 
new definitions are constructed to be directly operational from the experimental 
point of view. Finally, they show why without the thermodynamic limit or at 
phase-transitions, the systems do not have a single peaked distribution in phase 
space. 

In a closely related contribution. Dieter Gross makes the statement, 
that the microcanonical ensemble with Boltzmann's principle S' = /cg In 12 is the 
only proper basis to describe the equilibrium of a closed " small" system. Phase- 
transitions are linked to convex (upwards bending) intruders of the entropy, 
where the canonical ensemble defined by the Laplace transform to the intensive 
variables becomes multi-modal, non-local, and violates the basic conservation 
laws. The one-to-one mapping of the Legendre transform being lost. Gross claims 
that it is all possible without invoking the thermodynamic limit, extensivity, or 
concavity of the entropy. 



3 Thermodynamics 

3.1 Inequivalence of statistical ensembles 

Following the example exhibited long time ago by Hertel and Thirring , it is 
striking that these systems could lead to inequivalences between microcanonical 
, canonical or grand canonical ensembles. In this book, the first example is given 
by Barre et al Q who present the Blume-Emery-Griffiths (BEG) model which 
allows a deep understanding of the fundamental reason why this happens. They 
studied the spin- 1 BEG model both in the canonical and in the microcanonical 
ensemble. The canonical phase diagram exhibits a first order and a continuous 
transition lines which join at a tricritical point. It is shown that in the region 
where the canonical transition is first order, the microcanonical ensemble yields 
a phase diagram which differs from the canonical one. In particular it is found 
that the microcanonical phase diagram exhibits energy ranges with negative 
specific heat and temperature jumps at the transition energies. The global phase 
diagrams in the two ensembles and their multicritical behavior are calculated and 
compared. 

Pierre-Henri Chavanis shows similar features in self-gravitating systems 
where canonical and microcanonical tricritical points do not coincide either, as 
shown in Fig. ^ in the framework of self-gravitating fermions. Let us empha- 
size that this property survives to the introduction of a finite cut-off instead of 
quantum degeneracy as discussed by Chavanis. 
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Fig. 4. Inverse temperature as a function of the energy for self-gravitating fermions 
without cut-off. CE and ME refers respectively to the tricritical points in the canonical 
and microcanonical ensembles. 



3.2 Negative specific heats 

This fact produces striking phenomena in the microcanonical ensemble, since 
it may result in a negative specific heat region, as was emphasized by Ed- 
dington in 1926 and then discussed by Lynden-Bell |jl^. A first remark on the 
possibility of having negative specific heat in the microcanonical ensemble can 
even be found in the seminal paper on statistical mechanics by J.C. Maxwell ||l9|] . 
Thirring [p^ has finally clarified the point by showing that the paradox disap- 
pears if one realizes that only the microcanonical specific heat could be negative. 

Indeed, in the canonical ensemble the mean value of the energy of a system 
with different energy levels Ei is 



d\nZ 

8(3 



(16) 



where Z is the partition function. It is then straightforward to compute the 
specific heat 



C.^^^{{E-{E})'}>0 



(17) 



This clearly shows that the canonical specific heat is always positive. Notice also 
that this condition is true for systems of any size, regardless of whether a proper 
thermodynamic limit exists or not. 

This is not the case if the energy is constant as shows the simplified following 
derivation for the example of interacting self-gravitating systems. Using the virial 
theorem for such particles 



2{E,) + (Epot) = 



(18) 
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one gets that the total energy 

E = {E,) + {Epot) = -{Ec) ■ (19) 

As the kinetic energy E,, is by definition proportional of the temperature one 
gets that 

Loosing its energy, the system is becoming hotter. 

It is important at this stage to make a short comment on the Maxwell con- 
struction, usually taught in the framework of the Van der Waals liquid-gas tran- 
sition. The existence of a negative specific heat region corresponds to a convex 
intruder in the entropy-energy curve, as shown in Fig. ^. When the interactions 
are short range, the system will phase separate in two parts, corresponding to the 
two phases 1 and 2 with a molar fraction a;, so that the free energy xFi + {\—x)F2 
is lower than the original free energy. This is clearly possible if the energy cost 
of the interface is proportional to the surface whereas the energy gain is pro- 
portional to the volume of the phase. However, this is not any more possible 
when the interactions are long range since, on one hand, it is not straightfor- 
ward to define a phase and, moreover, the system is not additive. The Maxwell 
construction has to be redefined in this new framework. 




Fig. 5. Schematic shape of the entropy S as a function of the energy E with a convex 
intruder: the solid curve corresponds to the microcanonical result, whereas the dashed 
line to the canonical one. 



Let us note that the microcanonical entropy as a function of the energy and 
of the order parameter generically leads to the landscape presented in Fig. ^ 
The projection for the critical points of the surface onto the entropy-energy plane 
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Fig. 6. A stylized microcanonical entropy as a function of the energy and of the order 
parameter mimicks an Alpine landscape where the workshop took place. The projection 
for the critical points of the surface onto the entropy-energy plane produces the well 
known "swallowtail" catastrophe. 



produces the well known " swallowtail" catastrophe [gO| , depicted on the right of 
the figure. 

This concept of negative specific heat is now widely accepted in the as- 
trophysical community, and was popularized in particular by Hawking pl| | in 
1974, with some esoteric applications to black holes. The caloric curve of self- 
gravitating fermions derived by Chavanis and shown in Fig. ^ emphasizes such 
negative specific heat branch: the dotted branch is one example. Similarly one 
gets negative specific heat branch in the BEG model proposed by Barre et a/ Q . 
In the canonical ensemble, they correspond to local maxima or saddle point of 
the corresponding free energy; it is the constraint of keeping the energy constant 
that stabilizes these canonical unstable states in the microcanonical ensemble. 

Experimental groups have recently claimed signatures of negative specific 
heats in small systems. The first one corresponds to nuclear fragmentation 
even if the authors use prudently the word "indication" of negative specific heat. 
The latter being inferred from the event by event study of energy fluctuations 
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from Au + Au collisions. However, the signatures correspond to indirect mea- 
surements. 

In the clusters community, two experimental groups have very recently re- 
ported negative specific heat. The first system correspond to atomic sodium 
clusters, namely Na^4y and the negative microcanonical specific heat has been 
found near the solid to liquid transition. The cluster ion are produced in a gas 
aggregation source and then thermalized with Helium gas of controlled temper- 
ature. Accelerated thanks to the charge in a mass spectrometer, they are finally 
irradiated by a laser to determine the energy from the evaporation of several 
atoms after laser irradiation, also called photofragmentation. However, the con- 
trol of equilibrium is as always the key point and therefore the evaluation of the 
temperature seems to be questionable, in particular since the temperature could 
not be constant during the motion of the ions. 

In the Lyon's molecular cluster experiment |Q, with H^y, the energy and the 
temperature are determined from the size distribution of fragments after colli- 
sion of the cluster with a Helium projectile. To simplify the method, the larger 
the ratio of small fragments versus large ones, the larger is the temperature de- 
termined using the Bonasera et al procedure . The caloric curve reported Q 
shows a plateau. Work along this line is in progress and seems to show a negative 
specific heat region. 



3.3 Non extensive statistics 

Constantino Tsallis, Andrea Rapisarda, Vito Latora and Fulvio Baldovin |^ 
review the generalized non-extensive statistical mechanics formalism and its im- 
plications for different physical systems. The original very interesting idea is to 
generalize Boltzmann's entropy by defining 

S, - kB ^ ~ (21) 

q — 1 

where X^iPi = 1- Using either the L'Hopital rule or a first order expansion of 
the term pj in power of q, one immediately notices that 

lim Sq = -fcs y^Pi In Pi , (22) 
g— ^1 ^ — ^ 

i 

i.e. the well known Shannon entropy, known to be equivalent to the Boltzmann's 
one. 

However, for q different from 1, this generalized entropy Sq is non additive, 
and one gets 

Sq{A + B)^Sq{A)+Sq{B) + {l-q)Sq{A)Sq{B) . (23) 

They illustrate in particular its application and the meaning of the entropic in- 
dex q for conservative and dissipative low-dimensional maps. They also report 
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on non Boltzmann-Gibbs behavior p6[ and hindrance of relaxation for Haniiho- 
nian systems with long-range interaction, where fingerprints of the generalized 
statistics have recently emerged. 

This very interesting proposal |^ had however until now no strong founda- 
tions and many physicists were not ready to admit that the exponential Boltz- 
mann distribution of states is at equilibrium only a particular case of a gen- 
eralized distributions, with power tails. Dieter Gross in particular makes 
different comments to this point. On the contrary, Tsallis et al emphasize also 
different situations were the Boltmann-Gibbs behavior is clearly not appropriate. 

Recently, Beck and Cohen Q showed that considering different statistics 
with large fluctuations, one can obtain generalized results, called superstatistics, 
with the Tsallis formalism being presumably so far the most relevant example. 
Moreover, Baldovin and Robledo worked out |2^,^ exactly the q indices for the 
generalized largest Lyapunov exponent proposed by Tsallis. This an important 
step toward the derivation of a complete theory which, in particular, should help 
to understand the limits of its applications. 



4 Dynamical aspects 

An essential peculiarity of these physical systems, and of some of their simplified 
models, is that a classical system of particles with long range interactions will 
display strong non-equilibrium features. Dynamics is typically chaotic and self- 
consistent, since all particles give a contribution to the field acting on each of 
them: one calls this self-consistent chaos. Numerous physical systems fall in this 
category: galactic dynamics, dynamics of a plasma, vorticity dynamics,.... 

It is therefore essential to study the thermodynamic stability of these systems 
and in particular to understand the formation of structures trough instabilities. 
They should have logical similarities with the Jean's instability of self-gravitating 
systems, or with the modulational instability, leading to the formations of lo- 
calized structures, as confirmed by preliminary results. Additional dynamical 
effects, like anomalous diffusion and Levy walks, which are reported in the neg- 
ative specific heat regions, should be linked to these uncommon characteristics 
of thermodynamics ]30[| . 

In particular, Diego Del Castillo Negrete jo) discusses a mean-field single- 
wave model that describes the collective dynamics of marginally stable fluids 
and plasmas. He shows thus the role of self-consistent chaos in the formation 
and destruction of coherent structures, and presents a mechanism for violent 
relaxation of far from equilibrium initial conditions. The model bears many sim- 
ilarities with toy-models used in the study of systems with long range interactions 
in statistical mechanics, globally coupled oscillators, and gravitational systems. 

One of these toy models is for example studied by Dauxois et al [^. They 
consider the dynamics of the Hamiltonian Mean Field model which displays 
several interesting and new features. They show in particular the emergence of 
collective properties, i.e. the coherent (self-consistent) behavior of the particle 
motion. The space-time evolution of such coherent structures can sometimes be 
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understood using the tools of statistical meehanies , otherwise can be a manifes- 
tation of the solutions of an associated Vlasov equation. Both cases in which the 
interaction among the particles is attractive and the one where it is repulsive 
are interesting to study: they offer different views to the process of cluster for- 
mation and to the development of the collective motion on different time-scales. 
The clustering transition can be first or second order, in the usual thermody- 
namical sense. In the former case, ensemble inequivalence naturally arises close 
to the transition. The behavior of the Lyapunov spectrum is also commented 
and the 'universal' features of the scaling laws that it involves. 

Yves Elskens |^ shows that plasmas are a most common example of sys- 
tems with long-range interactions, where the interplay between collective (wave) 
and individual (particle) degrees of freedom is well known to be central. This 
interplay being essentially non-dissipative, its prototype is described by a self- 
consistent Hamiltonian, which provides clear and intuitive pictures of fundamen- 
tal processes such as the weak warm beam instability and Landau damping in 
their linear regimes. The description of the nonlinear regimes is more difficult. In 
the damping case, new insight is provided by a statistical mechanics approach, 
which identified the distinction between a trapping behavior and linear Landau 
behavior in terms of a phase transition. In the unstable case, the model has 
shown that the commutation of long-time and large- A'^ limits is not guaranteed. 

Chavanis considers also dynamical aspects in the framework of stellar systems 
and two-dimensional vortices. He discusses in particular two possible relaxation 
scenari: one due to collisions (or more precisely to discrete interactions) and the 
second one, called violent relaxation, really coUisionless but due to the mean 
field effect and the long range of the interaction. 

Finally, the dynamical processes that give rise to power-law distributions and 
fractal structures have been studied extensively in recent years. Ofer Biham and 
Ofer Malcai js^ describe recent studies of self-organized criticality in sandpile 
models as well as studies of multiplicative dynamics, giving rise to power-law 
distributions. Sandpile models turn out to exhibit universal behavior while in 
the multiplicative models the powers vary continuously as a function of the 
parameters. They consider the formation of a fractal object in the presence of 
a dynamical mechanism that generates a power-law distribution and presents a 
model that demonstrates clustering when the probability of adding a particle 
decays with a power a > d, so it has a short-range nature. 

5 Bose-Einstein condensation 

Finally, we would like to put a special emphasis on Bose Einstein Condensation 
(BEC), predicted by Bose and Einstein in 1924, which could be an important 
field of applications. With the recent achievement |Q of Bose-Einstein conden- 
sation in atomic gases thanks to the evaporation cooling technique, it becomes 
possible to study these phenomena in an extremely diluted fluid, thus helping to 
bridge the gap between theoretical studies, only tractable in dilute systems, and 
experiments. In the BEC, atoms are trapped at such low temperatures that they 
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tumble into the same quantum ground state creating an intriguing laboratory 
for testing our understanding of basic quantum phenomena. 

First, Jean Dalibard |Q presents how coherence and superfluidity are hall- 
mark properties of quantum fluids and encompass a whole class of fundamental 
phenomena. He reviews several experimental facts which reveal these two re- 
markable properties. Coherence appears in interference experiments, carried out 
either with a single condensate or with several condensates prepared indepen- 
dently. Superfluidity can be revealed by studying the response of the fluid to a 
rotating perturbation, which involves the nucleation of quantized vortices. 

Second, Ennio Arimondo and Oliver Morsch present the current investi- 
gations of Bose-Einstein condensates within optical lattices, where the long range 
interactions are an essential part of the condensate stability. Previous work with 
laser cooled atomic gases is also briefly discussed. 

On the theoretical side, the fluctuations of the number of particles in ideal 
Bose-Einstein condensates within the different statistical ensembles has shown 
interesting differences. Martin Holthaus explains why the usually taught 
grand canonical ensemble is inappropriate for determining the fluctuation of the 
ground-state occupation number of a partially condensed ideal Bose gas: it pre- 
dicts r.m.s. -fluctuations that are proportional to the total particle number at 
vanishing temperature. In contrast, both the canonical and the microcanoni- 
cal ensemble yields fluctuations that vanish properly for the temperature going 
toward zero. It turns out that the difference between canonical and microcanon- 
ical fluctuations can be understood in close analogy to the familiar difference 
between the heat capacities at constant pressure and at constant volume. The 
detailed analysis of ideal Bose-Einstein condensates turns out to be very helpful 
for understanding the occupation number statistics of weakly interacting con- 
densates. 

Ulf Leohnardt |Q shows that Bose-Einstein condensates can serve as lab- 
oratory systems for tabletop astrophysics. In particular, artificial black holes 
can be made (sonic or optical black holes). A black hole represents a quantum 
catastrophe where an initial catastrophic event, for example the collapse of the 
hole, triggers a continuous emission of quantum radiation (Hawking radiation). 
The contribution summarizes three classes of quantum catastrophes, two known 
ones (black hole, Schwinger's pair creation) and a third new class that can be 
generated with slow light. 

Finally, Gershon Kurizki presents an exciting theoretical idea to in- 
duce long-range attractions between atoms that acts across the whole Bose- 
Einstein condensate. He shows that dipole-dipole interatomic forces induced by 
off-resonant lasers 



^ {cos qr + qr SYR qr) q cos qr 



(24) 



allow controllable drastic modifications of cold atomic media. " Sacrifying strength 
for beauty" , Kurizki proposed |40j to average out the first term in 1 /r^ of the 
dipole-dipole interaction by the different lasers, in order to keep only the last one 
with a 1/r interaction. The important point is that induced gravity- like force 



Dynamics and thermodynamics of systems with long range interactions 



17 



would be strong enough to see it acting among atoms in the BEC: i.e. that, 
having induced the gravity-hke attraction in the BEC, one could switch off the 
trap used originally to create the BEC, and it will remain stable, holding itself 
together. 

Depending on the number of lasers, the resulting gravity- like force could be 
anisotropic for three lasers, or strictly identical to gravity with eighteen lasers ! If 
the last proposal is presumably too speculative and if the difficulties (the power 
of the laser required being really huge) facing the experimentalists are a real 
challenge, the ability to emulate gravitational interactions in the laboratory is of 
course fascinating. Indeed, these modifications may include the formation of self- 
gravitating "bosons stars" and their plasma-like oscillations, self-bound quasi- 
onedimensional Bose condensates and their "supersohd" density modulation, 
giant Cooper pairs and quasibound molecules in optical lattices and anomalous 
scattering spectra in systems of interacting Bosons or Fermions. These novel 
regimes set the arena for the exploration of exotic astrophysical and condensed 
-matter objects, by studying their atomic analogs in the laboratory. 

6 Conclusion 

The dynamics and thermodynamics of long range system is a rich and fascinating 
topic in particular for the following issues: 

• Statistical physics: inequivalence of ensembles, negative specific heat, coUi- 
sionless relaxation, role of coherent structures, relationship between dynam- 
ics and thermodynamics, nonadditivity, generalization of the entropy,... 

• This problem has also the nice property to be related to neighbooring sci- 
entific disciplines. Not only mathematics with application of catastrophe 
theory |3^, large deviations theory but also to computer science; be- 
cause of the long range interactions naive numerical code are of order 7V^, 
and need the developments of efficient algorithms such as the heap based 
procedure |Q, local simulation algorithm for Coulomb interaction |43[,... 

• However, this methodological and fundamental effort should provide a gen- 
eral approach to the problems arising in each specific domain which has 
motivated this study: astrophysical objects, plasmas, atomic and molecular 
clusters, fluid dynamics, fracture, Bose-Einstein condensation, ... in order to 
detect the depth and the origin of the observed analogies or, on the contrary, 
to emphasize their specificities. 

Many of these different aspects are considered in this book but it is clear 
that, rather than closing the topic, it opens the pandora box. 
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